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ABSTRACT 



Confexf. The recent finding of substantial masses of cold molecular gas as well as young stellar populations in the host galaxies of intermediate 
luminosity quasars is at odds with results of Hubble Space Telescope imaging studies since the latter appear to yield mature, quiescent early 
type hosts. 

Aims. The characterization as 'quiescent' for the host galaxies of QSOs is being addressed. 

Methods. The radio and far-infrared properties of both the HST sample and a larger comparison sample of uv-excess selected radio-quiet 
QSOs are being analyzed. 

Results. Consistency is found with substantial recent or ongoing star-formation activity. 



Key words, quasars: general - galaxies: starburst 



1. Introduction 

Quasi-stellar objects, or QSOs, represent episodes of extreme 
nuclear luminosity in galaxies, over a wide wavelength range. 
QSOs display a considerable dispersion in range-loudness, 
which is commonly expressed (Kellermann et al. 1989) as ra- 
tio R of cm-radio over blue (4400A) flux density. Radio-quiet 
QSOs typically have R < 1, while radio-loud QSOs are charac- 
terized with R in excess of 10. QSO host galaxies - both of the 
radio-loud and the radio-quiet subgroup - have been targets of 
studies for many years. Early investigations (e.g., Hutchings et 
al. 1984) revealed luminous hosts, frequently with close com- 
panions, and sometimes displaying spectroscopic features of 
young stars. Systematic studies of QSO host galaxies became 
possible with the refurbished Hubble Space Telescope (HST): 
host galaxies of different nature were identified (Bahcall et al. 
1997). Using bulge-black hole mass correlations, these studies 
permitted assessment of the QSO accretion rates, and rates up 
to 20 per cent of the Eddington values were inferred (McLeod 
et al. 1999). Ground-based near-infrared studies of samples 
of both radio-loud and radio-quiet hosts as well as powerful 
(Pradio > 10 27 W/Hz) radio galaxy hosts indicated (Taylor et 
al. 1996) the predominance of very luminous spheroidal hosts, 
particularly in the two radio-loud AGN populations. These re- 
sults were confirmed by imaging studies from space, using 
HST. In the final paper (Dunlop et al. 2003) in a series dealing 
with WFPC2 studies of the host galaxies of matched samples 
of radio-loud and radio-quiet QSOs as well as radio galaxies, 
it was concluded from 2-D model fitting procedures that "... 
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for nuclear luminosities brighter than V = 23.5 the hosts are 
virtually all massive ellipticals with properties indistinguish- 
able from those of quiescent, evolved low-redshift ellipticals of 
comparable mass." 

Powerful radio-galaxies are long known to be hosted by 
massive ellipticals, which are however sometimes of disturbed 
morphology (e.g., Heckman et al. 1986). While the HST imag- 
ing results (Dunlop et al. 2003) seem consistent with AGN 
unification models (involving obscuring dusty circumnuclear 
tori - e.g., Barthel 1989, Urry & Padovani 1993) it must be 
stressed that also the presence of extended dust, cold gas, 
and young stars in many radio-loud AGN has been reported 
(Hes et al. 1995, Martel et al. 1999, Tadhunter et al. 2005, 
Van Gorkom et al. 1989), which in turn is difficult to recon- 
cile with the picture of quiescent ellipticals. Similar contrast- 
ing views exist on the hosts of radio-quiet QSOs, i.e., QSOs 
displaying radio-optical flux ratio's of order unity or smaller. 
Whereas - as mentioned above - 2-D model fitting of the 
red starlight distributions generally seems to yield dominant 
spheroids, the properties of the interstellar media of several 
radio-quiet QSO hosts do not support the quiescent massive 
elliptical classification. The recent finding (Evans et al. 2001, 
Scoville et al. 2003) of substantial masses of cold molecular 
gas, as well as several cases (Canalizo & Stockton 2001) of 
young circumnuclear stars in radio-quiet QSO host galaxies is 
at odds with conclusions of the imaging studies. Some radio- 
quiet QSOs must be (Clements 2000) post-starburst systems, 
of which QSO UN J 1025^-40 provides a spectacular exam- 
ple (Brotherton et al. 1999). An extensive radio-optical study 
(Miller et al. 1993) using 89 low redshift QSOs, indicated that 
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roughly half of the weak radio emission in the radio-quiet sub- 
group might draw from host star-formation. A study of the in- 
tegrated radio and far-infrared emission in various classes of 
active objects (Colina & Perez-Olea 1995) confirmed that sus- 
picion. The HST results (Dunlop et al. 2003) are also inconsis- 
tent with modeling studies (Rowan-Robinson 1995) of infrared 
spectral energy distributions (SEDs), concluding "... an AGN 
is inevitably accompanied by a star-burst ..." 

We therefore examine in detail radio imaging and far- 
infrared photometric data for the radio-quiet subset of the HST 
QSO sample of Dunlop et al. (2003), and demonstrate from 
comparison with similar data for the uv-excess selected PG 
QSOs that the radio and far-infrared properties of both samples 
are in fact consistent with the presence of gas-rich star-forming 
interstellar media in the AGN hosts. So while we acknowledge 
the presence of a massive spheroid as a necessary condition 
for the nuclear activity to occur, we conclude that the circum- 
nuclear regions in the host galaxies of intermediate luminosity 
QSOs must contain young stars, supernova remnants, and star- 
burst heated dust. That conclusion is also in line with results 
of recent SDSS studies (Kaufmann et al. 2003), suggesting that 
many AGN are post-starburst systems, and recent ISO studies 
(Haas et al. 2003). 



2. Sample and analysis 

The HST/Dunlop et al. radio-quiet (i.e., radio-weak, not radio- 
silent) QSO sample (Dunlop et al. 2003) consists of thirteen 
objects, in the redshift range 0.1 < z < 0.25. Of these thirteen, 
QSO 1635+1 19 is associated with an extended (~1 arcmin) ra- 
dio source of luminosity P(5 GHz) = 5 x 10 24 W/Hz; combining 
this with its R(5 GHz) value of 67, we contend that this object 
should be formally classified as a radio-loud AGN (its optical 
magnitude is too faint for a QSO). Note that radio-emission in 
radio-quiet QSOs is generally weak (in absolute and relative 
terms), and confined to the nuclear and/or circumnuclear re- 
gions of the QSO host galaxy. The thirteen HST QSOs will be 
compared with Palomar-Green (PG - Schmidt & Green 1983) 
QSOs in the same redshift range, 0.1 < z < 0.25. The latter is a 
well known sample of QSOs selected through their ultraviolet- 
excess, and has well-documented properties (Sanders et al. 
1989): thirty-one radio-quiet PG QSOs are found in the rel- 
evant redshift range. Both samples are tabulated in the Table 
(note: six QSOs in common). Besides name, redshift, and ab- 
solute optical magnitude, we list flux densities at 60yum and 
5 GHz, inferred luminosities at 60/vm and 5 GHz, radio core 
fractions f c at 5 GHz, u-parameters (logarithmic ratio's of the 
60yum and 5 GHz flux densities), for the integrated (u) and 
the extended, radio-core subtracted (u') emission, and rele- 
vant references. Core fractions f c were obtained by compar- 
ing the 5 GHz flux densities obtained at low resolution (typi- 
cally 18 arcsec - Kellermann et al. 1989) with those at high 
resolution (< 0.5 arcsec). In several cases IRAS ADDSCAN 
routines were used to re-assess the 60yum properties, yield- 
ing new detections as well as improved upper limits. Having 
-25 < M fl /,5 < -23, the objects comprise a representative sam- 
ple of intermediate luminosity radio-quiet QSOs. 
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Fig. 1. 60yum FIR luminosities (expressed in log L ) versus ab- 
solute (nuclear) optical magnitudes, for the QSOs from both 
samples (with the exception of 0257+024, 1635+119 and 
2344+184, whose optical faintness formally disqualifies these 
object as QSOs). The crosses mark QSOs from the HST sam- 
ple, the open hexagons mark QSOs from the PG-sample. For 
the sake of simplicity, the R-values for the HST sample are 
combined with the B-values for the PG-sample (R-values used 
for the common objects), since the (small) color transforma- 
tions are uncertain due to intrinsic spread in the optical contin- 
uum slopes and varying emission line contributions. 



3. Results 

The 60/im far-infrared luminosities, for the HST and the PG 
samples, run from 10.5 to 12.5 (in log L Q ). The former has five 
upper limits, the latter eight. The radio luminosities run from 
10 22 to 10 24 W/Hz, with one 10 24 7 in the Dunlop set (the radio- 
loud object). Although the fractions of 60/mi non-detections 
(5/13 vs. 8/31) suggest otherwise, survival analysis (taking the 
upper limits into account) cannot exclude that the two samples 
were drawn from one and the same parent population. More 
sensitive infrared photometry is needed to re-address that issue. 
The FIR and radio luminosities however, indicate that the hosts 
cannot be quiescent early type galaxies. The Dunlop and PG 
QSO sample display 60/im luminosities which are orders of 
magnitude in excess of (IRAS/ISO) luminosities for massive 
inactive ellipticals, which are generally (Walsh et al. 1989) in 
the range 10 7 - 10 9 L . There would be an obvious explanation 
for this fact, if all far-infrared and radio emission drew from the 
nuclear activity, but this appears not to be the case as we will 
demonstrate in the following. 

Fig- n com P ares the optical nuclear luminosities and the 
60/L<m luminosities, for the QSOs in both samples. While the 
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two quantities are broadly correlated, a substantial scatter is 
seen. Fig. [2- dealing with a small sample in a fairly narrow 
redshift range - confirms results from the Haas et al. (2003) 
study of a large sample of PG QSOs studied with ISO. This 
study reports the presence of high infrared luminosities imply- 
ing large dust masses ranging from 10 6 to 10 8 5 M . However, 
in contrast to the near-IR (3-10pm) and the mid-IR (10^40//m) 
luminosity, being fairly well correlated with the QSO (blue) lu- 
minosity, the far-IR emission (40-150ji/m) is found to display 
a scatter of about a factor 10 in its correlation with that QSO 
luminosity. This suggests that FIR-emission surely draws from 
the active galactic nucleus but must in addition have a signif- 
icant contribution from stellar heating. As Haas et al. (2003) 
stress, the large amounts of dust in QSOs would require that the 
QSO phase was at least preceeded by a starburst. We investi- 
gate this further by addressing the radio continuum emission as 
this offers an unobscured probe of host galaxy star-formation. 

As demonstrated earlier (e.g., Wilson 1988), radio- 
emission in active galaxies can be separated into two com- 
ponents, one deriving from the nuclear activity and one from 
host galaxy star-formation. Subtracting the AGN-driven com- 
ponent from the integrated radio emission, an improved cor- 
relation with respect to the integrated long wavelength far- 
infrared emission (which draws from star-formation) is being 
found. The u-parameter, as defined above, is generally used 
for this purpose. This parameter relates the 60/mi FIR emis- 
sion to the 20cm or 6cm radio emission. Star-forming ob- 
jects are characterized (Condon & Broderick 1988) with av- 
erage (20cm) parameter 1120 = 2.0, and ±2<x-range 1.6-2.4, 
i.e., 1.4 GHz radio emission with strength between 0.4 and 2.5 
per cent of the 60yum FIR emission (1 per cent in the mean). 
Using the standard radio spectral index -0.7 for optically thin 
synchrotron emission, this translates to 6cm (5 GHz) parame- 
ters in the range 1.95 < Ug < 2.75 (mean 2.35). In fact, from 
radio imaging and FIR SED decomposition it is well known 
that the FIR emission combines an AGN- with a star-formation 
component. Objects with an appreciable active nucleus display 
surplus radio emission, hence u-parameters < 1. In Seyfert-2 
galaxy NGC 1068 for instance, roughly 80% of its 20cm radio 
emission draws from the 1.5 kpc triple radio source, hence is 
generated by the AGN. The remaining 20% occurs on larger 
scales, and finds its origin in host star-formation (Telesco et 
al. 1984). Prototypical ultraluminous starburst galaxy Arp220 
(Sanders & Mirabel 1996) has u 6 = log S(60yum)/S(5 GHz) = 
log (103.3Jy/0.204Jy) = 2.70 (hence negligible AGN compo- 
nent), star-forming QSO Mkn231 (Carilli et al. 1998) has u 6 
= log S(60//m)/S(5GHz) = log (35.4Jy/0.419Jy) = 1.93 butu fi 
= log S(60/im)/S(5GHz, extended) = log (35.4/0.173) = 2.31. 
In words, roughly 40% of the 5 GHz radio emission in Mk 23 1 
(and most of the long wavelength FIR emission) must draw 
from star-formation. 

We therefore use sensitive, high resolution radio imaging 
data (references in Table column 12) to assess - where possi- 
ble - the radio-core fractions (/ e ) and to transform the global, 
integrated u-parameters into AGN-subtracted u' -parameters. It 
must be noted that we are dealing here with radio-quiet ob- 
jects, i.e., objects without large scale radio lobes or jets, and 
that weak AGN radio emission - if any - is most likely concen- 



trated in faint subarcsec scale radio components (cores here- 
after). Inspection of the f c entries shows several blanks in the 
case of objects lacking adequate radio imaging, but on the 
other hand very few pure-core (f c = 1) objects: diffuse galac- 
tic scale radio emission must be occurring often. Given that 
no cases of f c > 1 are seen, the f c values cannot be attributed 
to core variability combined with different observing epochs: 
they must reflect the presence of weak, diffuse, arcsec-scale 
(host galaxy scale) radio emission which is resolved out at 
high angular resolution. Typical diffuse luminosities are of or- 
der 10 22 - 10 23 W/Hz, i.e., roughly one to ten times the diffuse 
radio emission in our Galaxy. Inspection of the Table more- 
over shows that the extended (or integrated) radio emission 
frequently correlates with the integrated 60//m emission ac- 
cording to the radio-FIR correlation: many cases of U6 and u 6 
in the range 2.0 - 2.8 are seen. The few cases of exception 
(PG0923+201, PG0947+396, PG 1114+445, PG 1115+407, 
PG1435-067) display suspiciously high FIR emission which is 
probably due to contamination in the extended infrared beam. 
A few cases of mildly excessive radio emission (1.5 < U6 < 2) 
can be also seen. 

As examples we consider the QSOs PG 0052+251, 
PHL909 (0054+144), and PG 1012+008. Detailed HST im- 
ages of the first two objects have been presented (Bahcall et 
al. 1996), establishing a spiral host for the former and an ellip- 
tical host for the latter. PG 0052+251 displays several HII re- 
gions in the disk of its host, and the HST imagery (Dunlop et al. 
2003) indeed yielded a bulge/disk composite. As for PHL 909, 
the 2D-analysis (McLure et al. 1999) improved on the earlier 
(Bahcall et al. 1996) HST picture, excluding a disc host and 
establishing an R = -23.62 elliptical host having an 8 kpc 
half-light radius. As seen from the Table, both objects have 
significant diffuse host galaxy radio emission, and their AGN- 
corrected Ug-parameters fall nicely on the FIR-radio correla- 
tion. In particular, the luminous FIR emission of PHL 909 as 
well as the fact that 50% of its 6cm radio emission must orig- 
inate in extra-nuclear regions are noteworthy. While there can- 
not be any doubt on a disk+bulge host for PG 005 1+251, also 
the host of PHL 909 cannot be a quiescent mature elliptical. 
We recall that the residual images (McLure et al. 1999, Dunlop 
et al. 2003) of PHL 909 display faint excess emission towards 
the West; additional excess emission may well be present un- 
der the QSO PSF PG 1012+008 resides in a host galaxy which 
was found (Bahcall et al. 1997) to be in the process of inter- 
action with a neighbouring galaxy. Roughly one quarter of the 
radio emission has a diffuse arcsec scale origin, scaling with 
the 60/vm emission cf. the radio-FIR correlation. 

In other words, the radio/FIR differences with quiescent 
ellipticals cannot simply be attributed to the AGN. A sub- 
stantial fraction of the radio and FIR emission in these in- 
termediate luminosity QSOs must draw from star-formation 
in their hosts. While in some cases this star-formation can 
be readily connected to residual emission in the HST im- 
agery (after subtraction of the spheroids), in other cases it 
must connect to small residuals under the PSF, i.e., to star- 
formation within the central 100 parsecs. Given that circum- 
nuclear star-formation is frequently observed in Seyfert galax- 
ies (e.g., Gonzalez Delgado et al. 2001), in connection with 
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the fact that Seyfert galaxies are generally considered to be the 
low-luminosity counterparts of QSOs, these results are not un- 
expected. We stress that we are not questioning the presence 
of massive bulges in the QSO hosts (and the probable connec- 
tion (Dunlop et al. 2003) of the bulge/disk ratio to the QSO 
luminosity), but we feel urged to point out the special nature 
of the host ISM. Back of the envelope calculations indicate 
SFR's of order 10 M per year. It should be noted however 
that more luminous QSOs (brighter than -25) probably do not 
have scaled-up SFR's (Haas et al. 2003), unless at early epochs 
(e.g., Cox et al. 2005). Our results are nevertheless consistent 
with the multicolour imaging obtained by Jahnke et al. (2004a) 
for low redshift QSOs, as well as for more distant AGN by 
Jahnke et al. (2004b) and Sanchez et al. (2004). It is tempt- 
ing to interpret this and other evidence for circumnuclear star- 
formation within the AGN accretion model of Miralda-Escude 
& Kollmeier (2005). 

Finally, the above mentioned typical QSO SFR, computed 
from FIR data, is incompatible with the value obtained using 
the [OII]/}3727 emission line strength of these QSOs. This fact 
- pointed out by Ho (2005) - calls for some level of extinction 
towards the [Oil] emission, which is not unexpected given the 
large dust masses suggesting dusty young star clumps as f.i. 
also inferred in the Antennae galaxies. 

4. Summary 

We conclude that a QSO episode is not seldom accompanied 
or preceded by an episode of star-formation in the gas-rich cir- 
cumnuclear regions of its host galaxy. We anticipate that radio- 
loud QSOs will show the same effects, but note that radio imag- 
ing separating AGN and star-formation driven radio emission 
in these objects is not going to be simple. Decomposition of 
the FIR SEDS will be more promising in this respect: results 
from currently as well as soon flying infrared space observato- 
ries are eagerly awaited. The evolving dust and stellar content 
of QSO host galaxies will remain an issue of great importance 
and interest (e.g., Barthel & Sanders 2006). 
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Table 1. The Table lists integrated flux densities, in mJy, at 60/mi FIR and 5 GHz (6cm) radio, radio core fractions f c (fraction of 
5 GHz radio emission in subarcsec core), luminosities at 60/mi (in solar units) and spectral power at 5 GHz (W/Hz). Furthermore, 
u 6 -parameters, addressing conformity to the FIR-radio correlation for the integrated (u) as well as the core-subtracted emission 
(u'). Parameter U6 is defined as log S(60//m)/S(5 GHz), whereas u' 6 = log S(60//m)/S(5GHz, extended); star-formation generally 
yields u^-parameters in the range 2.0 - 2.8; radio-loud AGN (with strong non-thermal core and lobe emission) have U6 <k 2, incl. 
negative values. Values for the Hubble parameter and deceleration parameter of 50 and 0.5 respectively, in a A = universe have 
been used to compute source intrinsic quantities. The 13 QSOs from the HST sample (Dunlop et al. 2003) are followed by the 31 
PG (Schmidt & Green 1983) QSOs: the former have absolute (nuclear) R-band magnitudes listed, the latter B-band, or V-band 
(PG 0906+484). References: BWE91 - Becker et al. 1991, D03 - Dunlop et al. 2003, K89 - Kellermann et al. 1989, Ku98 - 
Kukula et al. 1998, S89 - Sanders et al. 1989, T96 - Taylor et al. 1996 
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a. inferred from 20cm integrated emission minus flat 6cm core spectrum 

b. IRAS Faint Source Catalog 



c. NRAO VLA Sky Survey l http : //www . cv . nrao . edu/nvss/ I 

d. own high resolution VLA measurements 

e. IRAS ADDSCAN analysis 

f . NASA Extragalactic Database l http : / /nedwww . ipac . caltech . edu I 



